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We  calculate  the  far  infrared  absorption  of  B  conformation  poly(dA) 
poly(dT),  using  eigenvectors  from  a  lattice  dynamical  model  refined  with  the 
latest  inelastic  neutron  scattering  data.  Our  various  calculations  include  the 
isolated  helix,  the  addition  of  counterions,  and  the  addition  of  a  "spine  of 
hydration"  in  the  minor  groove.  We  find  agreement  with  experimental 
infrared  data  except  for  one  mode.  We  then  attempt  to  model  a  coupling 
between  two  helices  and  show  some  evidence  that  this  mode  is  interhelical.  -  ■ 
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’  Since  our  initial  prediction  of  vibrational  modes  in  DNA  in  the  1-100  cm"1  r*  4  /*/ 

£ 
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frequency  range,  1  several  experimental  groups  have  used  Brillouin,  low  fre- 
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quency  Raman,  and  infrared  absorption  measurements  to  study  this  region 
of  the  spectrum  Resonant  rather  than  dissipative  modes  have  been  observed 
and  in  one  case  the  lifetime  of  such  a  mode  has  been  determined.  Infrared 
absorption  has  been  carried  out  at  both  room  temperature  and  at  low  tem¬ 
perature  The  relatively  sharp  lines  available  from  this  measurement  have 
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indicated  a  temperature  dependence  in  agreement  with  our  calculations.  The 
good  experimental  data  has,  however,  raised  new  challenges  to  the  theory. 

Can  one  achieve  good  fits  and  assignments  to  observed  lines  in  this  region? 

This  question  implies  fitting  both  the  frequency  and  the  relative  absorption. 

In  addition  to  this  new  Raman  and  infrared  data,  inelastic  neutron 

scattering  data  is  also  now  available. 1  This  neutron  data  is  of  particular 

interest  to  us  as  it  has  allowed  refinement  of  our  model  of  the  long  range  non- 
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bonded  interactions  in  DNA  These  non-bonded  interactions  are  the  predom- 

r  ■ 

inant  interaction  at  lowest  frequencies;  that  is,  v<  10cm'1.  The  non-bonded 
interactions  are  unimportant  at  higher  frequencies  (i/>  200cm  _))  where  the 
bonded  interactions  predominate,  yet,  they  are  still  important  for  determining 
frequencies  in  the  range  of  interest  (10  <  v<  200cm1)  in  this  paper.  Our  ear- - — - 
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The  low  frequency  dielectric  constant  of  water  is  very  large  but  the  higher 
frequency  value  is  considerably  less.  Lee,  et.al.^  have  used  acoustic  attenua¬ 
tion  data  to  estimate  the  rotational  relaxation  time  of  water  molecules  near 
the  DNA  helix.  They  indicate  that  the  nearest  layer  of  water  can  rotate  only 
up  to  frequencies  «  4  GHz,  and  that  bulk  water  can  rotate  up  to  frequencies  « 
80  GHz,  the  latter  in  agreement  with  other  estimates  of  the  rotational  relaxa¬ 
tion  time  for  bulk  water. 

The  5-20  GHz  acoustic  data  is  still  in  the  frequency  range  where  the 

water  dipoles  can  rotate  and  effectively  sheild  electric  fields.  The  inelastic 

neutron  data  is  at  frequencies  of  hundreds  of  GHz  and,  therefore,  above  the 

region  of  large  water  dielectric  response.  We  were  able  to  use  the  neutron 

data  to  fit  a  non- bonded  interaction  model  that  should  be  appropriate  for  even 

higher  frequencies,  where  the  principal  element  in  that  model  is  the  electros- 
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tatic  interaction  shielded  by  a  model  of  dielectric  response.  Once  above  the 
water  rotational  relaxtion  frequency  the  dielectric  constant  of  water  should  be 
a  relatively  slowly  varying  parameter.  This  paper  will  then  use  the  non- 
bonded  model,  based  on  fitted  inelastic  neutron  scattering  data,  along  with  a 
bonded  force  constant  model,  developed  by  refining  higher  freqency  Raman 
and  infrared  data^  ^  to  calculate  the  frequencies  and  the  absorption  of  lines  in 
the  far  infrared  in  poly(dA)  poly(dT)  to  compare  with  the  observed  far 
infrared  absorption. 

As  demonstrated  by  Powell,  et.al  relatively  sharp  infrared  lines  can  be 
found  in  good  samples  of  DNA.  In  these  best  samples,  the  DNA  is  more 
ordered  than  in  the  poorer  samples  which  give  rise  to  much  broader  absorption 
peaks  In  these  more  ordered  DNA  samples,  one  would  expect  strong  interheli¬ 
cal  interactions.  Such  strong  interhelical  interactions  bring  into  question  the 


calculation^  approximation  of  a  single  isolated  helix.  Previous  studies  have 
centered  on  the  intrahelical  modes,  associated  with  motions  mostly  in  a  single 
helix.  Now,  in  addition,  one  has  strong  likelihood  of  interhelical  modes  -- 
modes  in  which  the  interactions  between  helices  are  dominant. 

For  properly  analyzing  a  system  having  both  intra-  and  inter-helical 
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modes,  one  must  use  an  entire  crystalline  unit  cell.  According  to  Lindsay 
such  a  cell  would  have  at  a  minimum,  two  to  three-ten  base  pair  sections  of 
helix,  plus  all  the  counterions  and  water  of  hydration.  The  secular  matrix 
would  be  roughly  30x30  larger  than  our  homopolymer  matrices.  Methods  to 
approximate  such  calculations  are  being  developed. 

In  this  paper  we  present  a  simple  set  of  calculations  which  are  aimed  at 
separating  out  the  intra-  and  inter-helical  modes.  We  consider  how  well  a 
simple  homopolymer  lattice  dynamics  calculation  can  predict  the  frequency 
and  relative  absorption  of  the  intrahelical  modes.  We  do  this  by  comparing 
the  isolated  helix  vibrational  modes  to  those  from  several  calculations  in  which 
the  helix  interacts  with  its  surroundings  in  increasing  detail.  These  calcula¬ 
tions  are  for  interactions  with  counterions,  interactions  with  water  molecules, 
and  then  the  combination  of  the  two.  We  find  the  method  works  fairly  well 
for  those  modes  we  believe  to  be  intrahelical  Furthermore,  we  develop  a 
crude  calculation  to  predict  the  frequencies  of  the  helix  interacting  with 
another  helix  in  its  surroundings  --  a  calculation  of  interhelical  modes.  We 
find  that  the  principal  feature  in  the  measured  infrared  spectrum  that  does  not 
correspond  to  our  intrahelical  calculations  corresponds  roughly  to  a  frequency 
predicted  by  our  crude  interhelical  calculation 
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Calculations  of  Vibrational  Modes 

We  begin  with  the  application  of  theories  of  solid  state  physics  to  DNA 
and  model  the  DNA  homopolymer  as  a  one-dimensional  infinite  lattice  (or 
periodic  arrangement  of  atoms)  whose  periodic  units  (or  unit  cells)  are  nucleo¬ 
tide  pairs  or  base  pairs.  The  entire  double  helical  homopolymer  is  generated 
by  taking  one  base  pair  and  repeating  it  along  the  helix  axis  in  the  manner  of 
translations  and  rotations  about  this  axis  as  if  the  unit  cell  were  on  a  screw. 
In  fact,  such  a  lattice  is  said  to  have  a  screw  symmetry  with  a  translation 
operation  of  p,  the  pitch  of  the  helix,  and  a  rotation  operation  of  V,  the  pitch 
angle.  Each  application  of  this  screw  operation  yields  the  next  unit  cell. 
Then,  according  to  Bloch’s  theorem,  for  an  infinite  lattice,  the  motion  of  the 
atoms  in  one  repeat  unit  is  the  same  as  that  in  another  except  that  the  two 
may  be  of  different  phase.  Therefore,  the  problem  of  studying  the  entire  lat¬ 
tice  is  reduced  to  the  study  of  one  unit  cell  with  all  possible  phase  relation¬ 
ships  between  succeeding  unit  cells. 

Given  the  atomic  coordinates,  masses,  charges,  and  force  constants,  we 
find  the  kinetic  and  potential  energy  of  the  lattice,  the  equations  of  motion, 
and  thus,  the  secular  equation  as  in  Eyster  and  Prohofsky.^’^  The  secular 
equation  is  derived  in  terms  of  the  mass- weighted  cartesian  displacement  coor¬ 
dinates  (MWC)  7,: 

( 1 ) 


m,  being  the  mass  of,  and  *,  the  displacement  of  the  atom.  Then  the  secu¬ 
lar  equation  is  given  as 


|  B*(0)FB(ff)  I  1  o 


(2) 


j  ’j>  v  /  ,*  /  «•  /  / 
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B'(0)FB(0)  should  actually  be  written  as  IB^lFBf®)],  +  [B*(0)FB(0)]ia  ,where  the 
first  term  is  for  valence  or  bonded  force  constants,  and  the  second  term  is  for 
non-bonded  interactons. 

In  (B^flJFBftf)], ,  F  is  the  diagonal  matrix  of  force  constants  of  the  internal 

coordinates,  5,.  B(0)  serves  to  transform  F  to  MWC  coordinates  and  also 

account  for  the  phase  difference  8  betweeen  unit  cells  n,  B(0)  =  £  Bn  R"  e'*9, 

» 

where  R  is  a  rotation  matrix.  The  elements  of  B  are  calculated  as  in  Chapter 

15  11 

4  of  Wilson,  et.al.  The  valence  force  constants  are  refined  by  Lu,  based 

on  spectral  data  above  400  cm1.  Thus,  the  rows  of  B  correspond  to  the 
numbering  of  internal  coordinates  and  the  columns  to  the  numbering  of  the 
x,y,z  components  of  the  MWC  coordinates. 

For  [B^tflFBfff)]^  ,  the  same  concepts  hold  for  F  and  B,  although  in  prac¬ 
tice  we  store  the  matrix  elements  somewhat  differently.  We  treat  the  non- 
bonded  interactions  between  the  two  atoms  of  different  unit  cells  as  a  stretch 
internal  coordinate  Since  all  of  these  interactions  are  stretches,  we  do  not 
number  the  internal  coordinates  but  instead  label  them  according  to  the  atoms 
involved. 

As  an  example,  we  compare  B  for  the  bonded  term  with  B  for  the  non- 
bonded  term.  For  the  bonded  term,  BXh  corresponds  to  the  third  internal 
coordinate,  53  and  the  y  component  of  the  MWC  coordinate  of  the  second 
atom  (atom  1  components  are  1,  2,  3,  atom  2  components  are  4,  5,  6,  ...).  If 

5  3  is  a  stretch  between  atoms  2  and  6,  B3A  ,  B  3  &  ,  B36  ,  Z?31f,  ,  i?317  ,  BJit  will  be 
the  only  non-zero  elements  on  row  three  But  in  the  non-bonded  term  B  is  a 
three-dimensional  matrix,  B  2  6  3  corresponding  to  the  y  (or  second)  component 
of  atom  2  (since  the  third  dimension  counts  x,y,z  for  the  atom  listed  in  the 
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first  dimension)  involved  in  the  stretch  with  atom  6.  Since  the  non-bonded 
problem  involves  around  (40)J  stretches  and  no  other  type  of  internal  coordi¬ 
nate,  the  three-dimensional  B  is  more  compact  for  computing  purposes. 


Upon  adding  these  parts  together,  we  can  solve  the  secular  equation  for 
the  eigenfrequencies  and  eigenvectors;  that  is,  the  resonant  frequencies  and 
the  relative  amplitudes  of  the  atoms  of  the  lattice.  Thus  we  have  an  adequate 
expression  of  its  dynamics. 

Modification  of  Non-Bonded  Potential  Energy 

The  potential  energy  comprises  bonded  and  non-bonded  terms  which 

leads  to  B’(0FB(0)  =  [B’(«)FB(«)1,  +  [B*(*)FB(0)]mt .  We  follow  the  method  of 
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Eyster  and  Prohofsky  in  calculating  the  non-bonded  force  constants,  modi¬ 
fying  the  functional  form. 

F,,  =  Ft  +  F,  (3) 


F  ,  =  2V 


/ - 

V  «,  (  ,  r 


3  ’ 
»  J 


(4) 


F  j  =  42 


(5) 


where  i  and  j  count  atoms  in  different  unit  (ells.  F,  is  for  long  range  or 
coulombic  interactions,  the  » •*  atom  having  a  charge  t,  ,  associated  with  a 
dielectric  screening  c,  and  separated  from  atom  j  by  a  distance  r, ,  The 
charges  are  taken  from  calculations  by  Miller.  the  hydrogen  atoms’  charge 


(as  well  as  mass)  being  included  in  the  atoms  to  which  they  are  attached.  The 
charges  are  reduced  by  a  factor  of  2.31  based  on  a  fit  to  inelastic  neutron 
scattering  data  as  well  as  comparison  to  a  few  charges  determined  experimen¬ 
tally  by  fitting  to  x-ray  scattering  data.^  r?=0.43,  which  includes  a  conversion 
factor  and  the  scaling  of  2.31  for  each  charge.  The  conversion  factor  is  needed 

because  the  units  of  the  force  constants  are  ,  but  the  charges  are  given 

in  units  of  electronic  charge  and  the  coordinates  in  angstroms.  Two  distinct 
dielectric  screenings  e,,«;  are  used,  as  some  atoms  (mostly  in  the  backbone) 
see  each  other  through  water,  whereas  others  (particularly  in  the  bases)  see 
each  other  through  DNA.  This  difference  explains  the  stronger  base-base 
interactions.  We  let  the  dielectric  vary  linearly  according  to  the  distance 
between  the  atoms  up  to  an  effective  separation  rtg  For  distances 

greater  than  rtg  the  dielectric  is  treated  as  a  constant.  If  the  atom  separation 
is  less  than  the  contact  distance  rcon  (the  minimum  separation  allowed 
between  the  centers  of  two  atoms),  the  dielectric  constant  is  taken  to  be  unity. 
We  use  an  effective  distance  of  10  A  and  a  contact  distance  of  2  A. 

The  formula  for  the  dielectric  function  is 


, 

+  ,  - - - rK; 

(  rtff  rron  ) 


(6) 


f  r  is  the  dielectric  constant  at  the  contact  distance,  rC(m  ,  ec  is  the  dielectric 
constant  at  the  effective  separation,  r,g  ,  and  r,;  is  the  separation  of  atoms  i 
and  j.  For  the  interaction  between  the  base  at"ms,  =  6  0  We  use  f„  =  9  0 


for  the  backbone  atoms’  interaction  As  discussed  in  the  introduction,  the 
choice  of  a  value  of  9  for  the  long  range  dielectric  value  of  water  near  DNA  is 


appropriate  to  calculations  in  the  THz  (30-300em  ~‘)  region.  F,  describes  the 
Van  der  Waals  interaction  when  n=8.  A  is  the  phenomonological  parameter 
adjusted  to  a  value  of  0.12  to  fit  the  infrared  absorption  data.  The  long  range 
interaction  is  for  the  unit  cell  and  nineteen  neighbors  on  either  side.  The  Van 
der  Waals  interaction  is  much  shorter  in  range,  involving  only  one  neighbor  on 
each  side  of  the  unit  cell.  The  new  set  of  eigenvectors  obtained  from  this 
refinement  are  used  to  calculate  the  absorption. 

Bare  DNA  Absorption  Calculation 

20 

We  apply  the  expressions  developed  by  Kohli,  et.al.  to  calculate  the 
absorption  of  radiation  by  poly(dA)  •  poly(dT)  which  is  in  the  B  conforma¬ 
tion.  For  the  relative  intensity  of  the  peaks,  we  calculate  the  Einstein  absorp¬ 
tion  coefficients: 

Bl  -  ^  I  <V«>  1 3  (7) 

bt  -  —  i  <*5>  r  (8) 

where  L  and  T  refer  to  the  longitudinal  (along  the  helix  axis)  and  transverse 
(perpendicular  to  the  helix  axis)  electromagnetic  fields  and  >  .>, 9>  is  the  transi 
tion  matrix  element  of  the  electric  dipole  moment  operator  for  the  1,9  vibra¬ 
tional  mode,  l  labeling  the  band  and  9  ,  the  phase  shift  from  cell  to  cell. 

Assuming  the  absorption  of  radiation  precipitates  a  transition  between 
the  two  lowest  vibrational  states,  these  matrix  elements  are  proportional  to 
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(10) 


the  sum  being  over  the  atoms  of  the  repeat  unit  with  charge  e,  and  mass  m, , 
with  N,  the  number  of  repeat  units  or  base  pairs.  Also, 


and 
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sin 


N  sin 


(ID 


(12) 


where  q‘t,  q,‘t,  are  the  components  of  the  /'*  vibrational  mode  eigenvector 
and  ^  is  the  angle  the  transverse  component  of  the  electromagnetic  field  makes 
with  the  f  direction 

For  our  case  of  an  infinite  chain,  the  absorption  occurs  only  for  9~  0.0  in 
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equation  9  and  9  tv  in  equation  10,  according  to  Kohli,  et.al  Then  for 
each  base  pair  (N- 'll. 
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Noting  that  w,tS_+^  =  w(|9__^  (due  to  time  reversal  symmetry),  and 

2, 

then  the  total  transverse  Einstein 

absorption  coefficient  is 
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Bt  =  (2) 
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(17) 


The  results  of  our  calculation  are  shown  in  figure  1.  We  have  plotted 
|  <n>  |  3  since  the  Einstein  coefficient  is  simply  a  multiple  of  that  value.  The 

5 

absorption  peaks  found  by  Powell,  et.al.,'  in  irradiating  Na  -  poly(dA) 
poIy(dT),  are  at  63,  83,  86,  100,  110  wavenumbers  (cm-1)  at  a  temperature  of 
7  K  as  shown  in  figure  2  Our  calculation  in  that  region  yields  sizable  peaks  at 
39,  45,  83,  89,  100,  109,  110  cm-1  for  9=  ip  and  100  cm-1  for  8=  0.0.  The  polari¬ 
zation  selection-,  are  different  for  the  9=  0.0  and  9=  ip  modes.  For  absorption 
at  8  0  0,  the  incident  electric  field  must  be  parallel  to  the  helix  axis.  For  the 
incident  field,  perpendicular  to  the  helix  axis  in  the  ordered  film,  the  absorp¬ 
tion  is  at  9=ip.  The  relative  contributions  to  each  depends  on  experimental 
details  We  see  no  peak  at  63  cm  1  and  do  not  have  a  mode  at  63  ,  m  1  How¬ 
ever,  outside  the  range  of  possible  observations  in  the  infrared  absorption 
experiment,  we  have  two  highly  absorbing  modes  at  39,  45  cm  1  1'sing  the 


12 


eigenvectors  of  the  modes,  we  find  the  largest  amplitude  motion  of  the  free 
phosphate  oxygens  occuring  in  the  39,  45  rm'1  modes  (the  motion  being  in  the 
adenine  backbone  for  39  cm-1  and  in  the  thymine  backbone  for  45  cm-1). 
These  modes  are  particularly  likely  to  be  affected  by  interhelical  interactions 
as  the  phosphate  groups  are  on  the  outside  of  the  helix  and  highly  charged.  In 
packing  two  right-handed  helices  parallel  to  each  other,  regions  of  close  con¬ 
tact  between  helix  backbones  necessarily  occur  as  they  must  cross  each  other 
or  distort  from  helical  symmetry.  It  is  likely  that  such  modes  are  pushed  up  in 
frequency  by  interhelical  interactions.  We  will  return  to  this  point  when  we 
discuss  our  calculation  with  interhelical  interactions. 

Counterions 

A  solution  of  DNA  requires  counterions  for  stabilization,  but  we  have  not 

included  the  presence  of  these  ions  explicitly  in  our  earlier  models.  Simula- 
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lions  by  Clementi  and  Corongiu  indicate  that  sodium  counterions,  strongly 
attracted  to  the  free  phosphate  oxygens,  follow  helical  patterns  along  the  back¬ 
bones  Such  interactions  with  the  phosphate  groups  could  alter  the  frequency 
of  phosphate  vibrational  modes. 

The  most  difficult  step  in  adding  counterions  to  our  model  is  determining 

tlu’ir  positions  relative  to  our  present  set  of  coordinates.  We  placed  sodium 

ions  an  aveiage  distance  of  2.7  A  away  from  the  free  phosphate  oxygens,  one 

counterion  for  each  phosphate  group.  This  distance  is  in  approximate  agree- 
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ment  with  calculations  by  Clementi  and  Corongiu.  Once  the  positions  are 
found,  t  lie  force  constants  can  be  calculated.  It  is  likely  that  the  main  interac¬ 
tion  of  the  ions  and  helix  is  via  electrostatic  forces.  Addition  of  counterions 
requires  an  increase  in  the  size  of  |B’( 0)FH( 9)}n,  .  the  non-bonded  term  in  the 


secular  equation.  The  new  elements  for  B  and  F  are  calculated  as  before  in 
the  non-bonded  case,  the  ions  interacting  with  the  DNA  atoms  sis  well  as  other 
ions  over  four  turns  of  the  helix. 

Solving  this  secular  equation  at  yields  eigenfrequencies  of  34,  40,  47, 

54  cm-1  in  the  range  of  interest,  34-65  cm-1,  as  compared  with  39,  45,  51,  54 

cm  1  for  the  bare  helix.  The  addition  of  counterions  yields  no  mode  near  60 

cm"1.  Calculation  of  the  absorption  with  these  eigenvectors  does  show  a 
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decrease  in  the  intensity  of  the  peaks  34-50  cm  l,  while  that  of  the  mode  at 
54  cm"1  and  the  modes  in  the  80  cm'1  region  showed  sharp  increases  in  absorp¬ 
tion.  The  exception  is  the  decrease  in  absorption  of  the  mode  at  74  cm'1.  In 
particular,  the  absorption  of  the  mode  at  89  cm'1  has  more  than  tripled,  while 
the  mode  at  83  cm'1  shows  only  a  slight  increase  in  absorption.  Instead  of  one 
peak  at  100  cm"1,  we  find  three  --  98,  99,  101  cm-1  -  and  none  at  110  cm'1. 
The  values  of  |  <p>  | 2  for  the  modes  in  the  range  of  50-115  cm'1  are  shown  in 
figure  3. 

Spine  of  Hydration 

The  water  molecules  of  the  solution  which  surround  the  atoms  of  DNA 

form  hydrogen  bonds  with  some  of  the  nitrogen  and  oxygen  atoms,  the  bond 

3 

lifetime  being  about  10  11  seconds  Thus,  for  frequencies  greater  than  a  few 
wavenumbers,  this  layer  of  water  molecules  is  part  of  (he  DNA  structure. 
Kopka,  et.al.  '  describe  the  most  distinct  pattern  of  water  molecules  in  B 
poly(dA)  poly(dT)  as  a  "spine  of  hydration".  In  her  data  the  spine  is 
formed  by  water  molecules  which  zig-zag  from  base  pair  to  base  pair  along  the 
I i*d i x  axis  in  the  minor  groove.  Between  two  base  pairs,  a  water  molecule 
forms  a  hydrogen-bond  with  a  thymine  02  atom  in  one  base  pair  and  an 


adenine  N3  atom  in  the  neighboring  base  pair.  All  along  the  helix,  water 
molecules  are  found  joining  neighboring  A-T  base  pairs.  Then  to  complete  the 
spine,  these  water  molecules  are  connected  to  each  other  by  forming  hydrogen 
bonds  with  another  water  molecule  located  between  those  which  bond  to  the 
bases  (figure  7  of  Kopka,  et.al.^). 

Aside  from  this  distinct  pattern,  a  total  of  about  20  water  molecules  per 
base  pair  are  found.  We  limit  ourselves  to  including  only  the  "spine  of  hydra¬ 
tion"  in  our  model,  the  difficulty  being  the  placing  of  the  molecules  in  our  unit 

cell  configuration  which  is  not  exactly  the  same  as  the  crystal  configuration  of 
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Kopka  et.al.  Then  we  calculate  the  modes  for  a  system  of  DNA,  ions,  and 
this  chain  of  water  molecules  interconnected  by  hydrogen  bonds  which  approx¬ 
imates  the  spine.  The  hydrogen  bond  lengths  for  our  calculation  vary  from 
2  64  to  2.89  A.  We  use  the  Lippincott-Schroeder  model^  to  calculate  the 
hydrogen  bond  stretch  force  constants  which  vary  from  0.66  (for  the  2.64  A 

bond)  to  0.08  (2.89  A  bond).  We  then  approximate  the  angle  bend 

force  constants  by  taking  the  average  of  the  stretch  force  constants  of  the 
bonds  involved  in  the  angle  bend  and  then  dividing  by  seven.  The  choice  of 

-  of  the  magnitude  of  the  stretch  force  constants  is  based  on  our  comparison 

of  the  valence  angle  bend  and  stretch  force  constants.  These  approximated 

angle  bend  force  constants  vary  in  strength  from  0.05  to  0.09  —  . 

rad  7 

Because  of  the  hydrogen  bonding  of  the  water  molecules  to  certain  DNA 
atoms,  we  increase  the  size  of  the  bonded  term  in  the  secular  equation, 
!B'(0)FB(0)), .  The  hydrogen  bond  force  constants  are  placed  along  the  diago¬ 
nal  uf  F  and  the  elements  of  B  are  calculated  as  before.  Including  the  mass  of 
the  hydrogen  atom  of  the  water  with  its  oxygen  atom,  we  carry  out  the  normal 


mode  calculation.  Using  these  new  eigenvectors,  we  then  calculate  the  absorp¬ 
tion.  Although  we  now  see  a  mode  in  the  60  cm-1  range  (at  67  cm-1),  it  shows 
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relatively  little  absorption.  In  the  range  70-110  cm-1  we  have  significant 
peaks  at  82,  86,  98,  110.  In  general,  the  absorption  has  decreased  relative  to 
the  bare  helix  calculation,  although  the  86  cm-1  mode  is  as  absorbing  as  the  89 
cm-1  mode  of  the  bare  helix  and  the  mode  at  110  cm-1  shows  an  increased 
absorption. 

Ions  Plus  Spine  of  Hydration 

Each  time  an  atom  is  added  to  the  model,  the  size  of  the  dynamical 
matrix,  Bf(0)FB(0)  increases.  Its  rows  and  columns  correspond  to  the  number¬ 
ing  of  the  MWC  coordinate  components  of  the  atoms.  Adding  ions  increased 
B'(0)FB(0)  through  the  non-bonded  term,  while  adding  water  molecules 
increased  B'(0)FB(0)  through  the  bonded  term.  Thus,  when  we  find  the  secular 
equation  for  this  system  by  adding  the  two  terms,  we  must  make  sure  the  ion 
and  water  matrix  elements  do  not  overlap.  The  addition  of  two  ions  and  two 
water  molecules  per  unit  cell  increases  the  size  of  the  final  matrix  for  the  unit 
cell  from  123x123  to  135x135.  The  eigenvectors  from  this  calculation  yield  a 

further  reduction  in  absorption  of  the  66  cm-1  mode  (from  0.25xl0-JO  to 
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0.03xl0-3°  ( e  -cm)3,  where  e  is  the  electronic  charge).  As  would  be  expected, 
in  general,  the  absorption  has  decreased  relative  to  the  calculation  in  which 
ions  alone  are  added  and  increased  relative  to  the  calculation  in  which  water 
molecules  alone  are  added.  Specifically,  in  the  70-110  cm  1  region,  the  modes 
with  significant  peaks  are  at  81,  86,  90,  97,  99,  101,  103  cm  1  as  shown  in 
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Although  the  addition  of  ions  does  seem  to  affect  the  modes  with  large 
amplitude  motion  of  the  free  phosphate  oxygens,  and  the  addition  of  water 
molecules  does  give  a  mode  near  63  cm-1,  such  results  are  still  not  in  good 
agreement  with  the  experimental  observation  of  a  highly  absorbing  peak  at  63 
fm'1.  More  likely,  some  other  mechanism  will  move  the  39,  45  cm-1  modes  to 
higher  frequency  --  namely,  interhelical  coupling. 

Interhelical  Coupling 

Another  mechanism  for  stiffening  of  the  39,  45  cm-1  modes  is  an  electros¬ 
tatic  interaction  between  the  phosphate  groups  of  different  helices  with  the 
sodium  counterions  included.  One  can  make  a  crude  attempt  to  model  such  a 
coupling  by  using  transverse  periodic  boundary  conditions.  We  essentially 
generate  an  infinite  array  of  base  pairs  (including  a  sodium  ion  at  each  phos¬ 
phate  group)  along  y  in  one  base  pair’s  (or  unit  cell’s)  frame  of  reference, 
where  y  is  approximately  parallel  to  the  interbase  hydrogen  bonds.  Then  we 
take  this  array  and  repeat  it  along  the  original  base  pair’s  helix  axis  by  the 
translations  and  rotations  that  give  DNA  its  helical  symmetry.  Of  course,  we 
have  not  succeeded  in  generating  separate  helices,  but  this  model  should  still 
be  valid  to  look  at  the  effect  of  the  interaction  of  phosphate  groups  of  nucleo¬ 
tides  of  different  helices  on  the  motion  of  the  free  phosphate  oxygens. 
Lindsay  provides  the  distance  between  helix  axes  from  fiber  diffraction 
experiments.  We  use  28  A  as  this  distance;  that  is,  in  generating  the 
transverse  array  of  base  pairs,  the  distance  between  the  helix  axis  of  each  base 
pair  is  28  A. 

We  use  the  same  electrostatic  formalism  as  for  the  calculation  of  the  addi¬ 
tion  of  counterions  to  the  model.  This  time  the  interaction  is  between  the  ±1 


transverse  repeat  unit  and  the  initial  repeat  unit  (cell  0),  where  ±1  refers  to  a 
translation  of  ±28Ay  .  Another  modification  is  due  to  the  amount  of  phos¬ 
phate  clashing  between  helices.  Our  model  includes  a  clash  for  each  phosphate 
group,  but  such  clashes  should  only  occur  approximately  every  five  vertical 
unit  cells  for  neighboring  helices.  For  this  reason,  we  decrease  the  force  con¬ 
stants  by  —  to  improve  the  frequencies,  realizing  that  details  of  our  calcula- 
tion,  such  as  amplitudes  of  specific  atomic  motion  will  be  less  meaningful. 

In  the  range  of  interest,  interhelical  coupling  yields  modes  at  33,  40,  45, 
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50,  52,  60,  61  cm  l.  These  modes  show  less  motion  for  the  free  phosphate 
oxygen  atoms.  The  modes  at  40,  45  cm-1  are  now  less  absorbing  than  the  39, 
45  cm-1  modes  of  the  bare  helix.  Also,  we  now  have  modes  with  higher 
absorption  in  the  50  -  65  cm'1  range  (figure  5),  whereas  before,  with  the  bare 
helix,  we  had  no  greatly  absorbing  modes  in  that  region.  With  the  addition  of 
ions,  only  one  significantly  absorbing  mode  appeared  in  that  range.  In  fact, 
the  absorption  of  the  52  cm'1  mode  is  near  the  magnitude  of  that  of  the  40 
and  45  cm'1  modes.  In  the  70-100  cm'1  range,  the  effect  on  the  peaks  relative 
to  the  bare  helix  is  similar  to  that  of  adding  counterions  alone,  except  that  the 
modes  are  shifted  slightly  higher:  84,  89,  98,  100,  103  cm'1. 

Conclusions 

We  do  find  that  our  single  helix  calculations  fit  fairly  well  the  observed 
modes  believed  to  be  intrahelical.  This  result  is  likely  as  the  modes  calculated 
in  this  frequency  regime  are  principally  stretch  modes  of  the  interbase  hydro¬ 
gen  bonds.  The  fact  that  even  the  addition  of  counterions  and  water 
molecules  shifts  these  modes  by  only  a  few  wavenumbers  reinforces  their 


effect  is  expected  to  be  particularly  important  in  those  cases  where  the  helices 
are  relatively  well  ordered.  Most  importantly,  these  well  ordered  cases  are  the 
ones  that  lead  to  the  best  resolved  spectra.  The  evidence  is  that  one  needs  to 
take  into  account  the  interhelical  interactions  to  explain  the  best  spectral  data. 
We  see  that  shifts  in  frequency  of  the  order  needed  to  obtain  better  agreement 
with  experiment  can  be  obtained  with  crude  approximations  of  the  interhelical 
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Figure  Captions 

Figure  1.  The  dotted  lines  represent  the  square  of  the  dipole  moment  matrix 
element  for  each  longitudinal  (9=  0.0)  mode  of  poly(dA)  •  poly(dT)  in  the  50 


-  115  cm-1  range.  The  solid  lines  are  for  the  square  of  the  dipole  moment 
matrix  element  for  each  transverse  (0=V>)  mode  of  poly(dA)  poly(dT)  in  the 
50  -  115  cm'1  range.  The  square  dipole  moment  unit  is  (e  cm  )5. 

Figure  2.  The  locations  of  the  peaks  in  the  measured  absorption  spectrum  of 
Powell,  et.al  ^ 

Figure  3.  The  absorption  is  given  as  in  figure  1  for  the  transverse  modes,  but 
with  the  addition  of  counterions. 

Figure  4.  The  absorption  is  given  as  in  figure  1  for  transverse  modes,  but  with 
the  addition  of  counterions  and  the  addition  of  water  molecules  binding  to 
atoms  in  the  bases  and  to  each  other 


Figure  5  The  absorption  is  shown  as  in  figure  1  for  transverse  modes,  but 
with  interhelical  coupling 
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